Stochastic kinetics of a single headed motor protein: 
dwell time distribution of KIF1A 
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KIF1A, a processive single headed kinesin superfamily motor, hydrolyzes Adenosine triphosphate 
(ATP) to move along a filamentous track called microtubule. The stochastic movement of KIF1A 
on the track is characterized by an alternating sequence of pause and translocation. The sum of 
the durations of pause and the following translocation defines the dwell time. Using the NOSC 
model (Nishinari et. al. PRL, 95, 118101 (2005)) of individual KIF1A, we systematically derive an 
analytical expression for the dwell time distribution. More detailed information is contained in the 
probability densities of the "conditional dwell times" t±± in between two consecutive steps each of 
which could be forward (+) or backward (-). We calculate the probability densities S±± of these four 
conditional dwell times. However, for the convenience of comparison with experimental data, we also 
present the two distributions of the times of dwell before a forward (+) and a backward (-) step. 
In principle, our theoretical prediction can be tested by carrying out single- molecule experiments 
with adequate spatio-temporal resolution. 

PACS numbers: 87.16.ad, 87.16.Nn, 87.10.Mn 
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I. INTRODUCTION 

Molecular motors are nano-devices which perform me- 
chanical work by converting part of the input energy; for 
the motors of our interest in this paper, the input is de- 
rived from the hydrolysis of ATP molecules |l] . In reality, 
these motors are also enzymes that hydrolyze ATP and 
utilize the input chemical energy to perform mechanical 
work. In this paper we specifically consider the mem- 
bers of a particular superfamily of motors, called Kinesin, 
which are involved in intracellular transport processes in 
living cells. This family is designated as KIF1A and the 
members of this family move along filamentous tracks 
called microtubule (MT) 0. 

One unique feature of a KIF1A is that, at least under 
the conditions of in-vitro experiments, it functions as a 
single- headed motor. The average properties, e.g., the 
average velocity, of these motors have been calculated 
analytically by using a theoretical model developed by 
Nishinari, Okada, Schadschneider and Chowdhury (from 
now onwards, referred to as the NOSC model) [3, H| 
which is an extension of the general approach pioneered 
by Fisher and Kolomeisky 

In single molecule experiments, individual motor pro- 
teins are observed to move in an alternating sequence of 
pause and translocation. The sum of the pause at a bind- 
ing site and the subsequent translocation can be defined 
as the corresponding "dwell time". Because of the in- 
trinsic irreversibility of mechano-chemical kinetics of the 
system, the inverse of the mean dwell time is the aver- 
age velocity of a motor. The dwell time distribution g(t) 
contains more detailed information on the stochastic ki- 
netics of a motor than that revealed its average velocity. 
For example, the randomness parameter 



< r > - < t > J 

<T>2 



(i) 



provides an estimate of the lower bound on the number 



of rate-limiting kinetic steps in each cycle of the motor 
6]. For some other motors, which move on nucleic acid 
strand, the analytical forms of the distributions of the 
dwell times have been reported recently @, HI ■ 

In this paper we report the exact analytical expression 
for the distribution of the dwell times of a KIF1A mo- 
tor in the NOSC model during a single processive run in 
between its attachment to the track and the next detach- 
ment. What makes the calculation more difficult in the 
case of KIF1A, compared to those of those reported in 
ref. 0, H| , is the occurrence of branched pathways in its 
mechano-chemical cycle. For motors which can step both 
forward and backward, one can define conditional dwell 
times which may be more easily extracted from the data 
obtained from single molecule experiments [3, [l(| . There- 
fore, in this paper we also report analytical expressions 
for the probability densities of these conditional dwell 
times as well as that of a few other closely related ran- 
dom variables. 



II. MODEL AND STOCHASTIC KINETICS 

A. The Model 

A MT track of the motor is modelled as a one- 
dimensional finite lattice having L number of discrete 
sites. Each site corresponds to a KIF1A binding site 
on the MT and the lattice spacing is the separation be- 
tween the successive binding sites on a MT. A KIF1A 
motor is represented by a particle with two possible 
chemical states labeled by the indices 1 and 2. The 
states 1 and 2 correspond to the strongly bound and 
weakly bound states, respectively, fig. [1] illustrates the 
detailed mechano-chemical cycle of KIF1A. The transi- 
tions lj ■<-> 2j are purely chemical whereas the transi- 
tions 2j O 2j±i, which correspond to the Brownian mo- 
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FIG. 1: Two state model for KIF1A. The indices ...,j — l,j,j + 
1, ... label the equispaced sites for the binding of the motor 
to its track. The states 1 and 2 correspond to the "chemical" 
states in which the motor is bound strongly and weakly, re- 
spectively, to the microtubule track. The allowed transitions 
are shown by the arrows along with the corresponding rate 
constants (transition probability per unit time). 



The Laplace transforms of the Fourier transforms 
of S(q,t) and W{q,t) are given by S(q,s) — 
r o °° dte~ st S(q,t) and W(q,s) = J °° dte~ st W {q,t). 
Therefore, using S(q,s) and W(q,s) as the two compo- 
nents of a vector P(g, s), the matrix equation (j4]) can now 
be written in a more compact notation as 

P{q, s) = R(g, s) _1 P(0), withR((j, s) = sl- M(q) (6) 

and P(0) is a vector whose elements are determined by 
the initial conditions for S and W. The initial conditions 
are S(j, 0) = and W(j, 0) = where k is an arbitrarily 
selected site. The dwell time at the site k is the total 
duration for which the motor stays at that site, starting 
from the initial condition mentioned above, irrespective 
of the chemical state, before its next departure from the 
same site. 

Thus the determinant of R(<?, s) is a 2nd order poly- 
nomial in s; that is [Til . 



|R|(?,«) = s 2 +a(q)s + -f(q) 



(7) 



where, 



tion, are purely mechanical. In contrast, the transition 
2j — > lj+i is mechano-chemical (see ref. [3] for a more 
detailed description). 



B. Kinetics and the master equations 

We define S(j,t) and W(j, t) as the probabilities of 
finding KIF1A in state 1 and state 2 respectively, at site 
j at time t. The master equations for these probabilities 
are 



dS(j,t) 
dt 



-u h S(j, t) + oj f W{j - 1, t) + u.W(j, t) (2) 



a(q) = uj h +2u b + u s + u f - uj b (e lq + e tq ) (8) 



7 (g) =w h (2u b +uj f )~u b uj h (e lq + e lq ) - u h u f e lq (9) 

Note that "f{q) is the determinant of the transition matrix 
M(q). Because of conservation of probability, in the q — > 
limit, all of the columns of M sum to zero. Therefore, 
we obtain |M(0)| = 7(0) = 0. _We define [Tl| the position 
probability density P(q, t) = S(q, t) + W(q, t), and hence 



s + a{0) 



■a(q)s + j(q) 



(10) 



C. Velocity and Diffusion constant 



dW(j,t) 
dt 



= w h S(j, t) -{u f + u) s + 2u b )W(j, t) 
+ u b (W(j-l,t) + W(j + l,t)) 



(3) 



with Wh — uj®[ATP], where [ATP] is the concentra- 
tion of ATP. Now we introduce Fourier transforms of 
S(j,t) and W(j,t) by S(q,t) = 

W(q,t) = YZL-oo W(j,t)e-™, where 1 = 
lattice spacing, d = 1. Thus, 



S(j,t)e- l <" and 
-1 and the 



S(q,t) 
W(q,t) 



M(q) 



S(q,t) 
W(q,t) 



(4) 



where the transition matrix M(q) is given by, 



M(q) 



(uj f e' tq +uj s ) 



uj h -{(2uj b + uj s +u f ) -uj b (e- iq + e lq )} 



(5) 



Following ref.[Tl[, the average velocity v of KIF1A and 
the diffusion constant D are found to be 



,7(0) = 

'a(O) (ujj + Lu h + lu s )' 



(11) 



and 



D 



7(0) - 2iva(Q) - 2f3{Q)v 2 _ 2u b oj h + w f u h - 2v 2 



2cv(0) 



2(u>f + uj h + lj s ) 



where 7(0) and 7(0) are the first and second derivatives, 
respectively, of 7(g) with respect to q evaluated at q = 
while /3(0) is the coefficient of s 2 at q — 0; in this 
case /3(0) = 1. Interestingly, the velocity v depends only 
on the coefficients 7 and a, of the lowest two orders of 
the polynomial obtained from the determinant of R(q, s). 
On the other hand, the diffusion coefficient D depends 
on the three lowest order coefficients /3, a, and 7 of the 
determinant of R(g, s). 
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FIG. 2: Dwell time distribution is plotted from equation (|30 
for few different values of W&. 



III. DWELL TIME DISTRIBUTION 

We define the individual forward and backward 
branching probabilities as p+,p and p- and the corre- 
sponding dwell time distributions are and ip-, 
respectively. Thus 



P+ 



L0 b 



2u b + Uf 



and, 



P- 



2u b + Uf ' 



u b 
2u b + Uf 



(13) 



(14) 



(15) 



We define the Fourier weights p+(q) = er %q and p-(q) = 
e lq for the forward and backward steps, respectively (step 
size d — 1 in our units) . Hence II lH , 



1 



sP(q,s) '{"m*' («)=»} 



1 - C(s)p+(q) - p-ip-(s)p-(q) 
l-C(*)-p-V-(«) 

(16) 

where £(s) — (p+^+( s ) + P ip (s)) and the symbol 
{pfc/fe' (?) = o} expresses the condition that the Fourier 

weight for all possible steps, except k is zero; in our case, 
k can represent either the forward (+) or the backward 
(— ) steps. For forward branching we obtain, 



sP(q,s) 



where 



«o 



\{p-(q)=o} 



1 - P+iP+(s)p+(<l) -Vi> (s)p+{g) 

1 -P+^+{s) -p'tp'(s) 
a + a+(s)p + (q) (17) 



s 2 + (2u b + Uh + Uf + u s )s + u h (2u b + Uf) 



s 2 + (uh + Uf + u s )s 



U b S + U h (u b + Uf) 

s 2 + (u h + Uf + u s )s 



(18) 



Thus, 



p+4>+{s) + pip (s) 



a 



U b S + U h (u b + Uf) 



s 2 + (2u b + u h + Uf + u s )s + u h (2u b + Uf) 
Inverse Laplace transform yields, 



(19) 



p+^+{t)+p^ (t) = 
where, 



Cie -nt/2 _ C2e -r*t/2 
2^ 



(20) 



r = J (2u b + u f + u h + u s ) 2 - A(2u b + u f )u h (21) 

Ci = {2ul + u b u f - u b u h - 2u f u h + u b u s + u b r Q ) (22) 

c 2 = {2ul + u b u f - u b u h - 2ufU h + u b u s - u b r ) (23) 

r% = 2u b + Uf + Uh + u s + r (24) 

r 2 = 2u b +Uf + Uh+u s — r (25) 



Similarly, 



1 



-|{ P+ ( 9 )=o} = a + a-(s)/5_(g) (26) 



sP(q,s) 

where ao is given by eq. (TT51) and 

UbS + u b u h 



a-(s) = - 



s 2 + {uf + Uh + u s )s 



(27) 



P-^-(s) = - 



a-(s) 



ao 

UbS + u b u h 



s 2 + (uh + 2ub + u s + Uf)s + Uh(2ub + Uf) 



(28) 



Inverse Laplace transform yields, 

u b {(n - 2u h )e— t ' 2 - (r a - 2u h )e~^l 2 } 



P-i>-(t) 



2r 



(29) 

Now total dwell time distribution can be written as fol- 
lows: 



g(t) = p + iP+(t)+ P Tp (t)+p_^_(t) 

{u b (r 2 - 2u h ) - UfU h }{e-^ 2 - e~^/ 2 ) 



ro 



(30) 



The total dwell distribution is plotted in the fig. [2j We 
observe that with the increase of Ub the most probable of 
dwell time shifts towards a smaller value. This is because 
rate constant Ub is related to diffusion of KIF1A in state 
2. For smaller values of Ub dwell time is dominated by 
the rate Uf, but with the increase of rate u b the dwell 
time depends on both u b and u f and eventually the mean 
dwell time decreases. 
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IV. PROBABILITIES OF SPLITTING AND 
CONDITIONAL DWELL TIMES 

A typical trajectory of a KIF1A motor consists of a 
random sequence of forward and backward steps. There- 
fore, we can define four different conditional dwell times 
t±±. Here t++ is the dwell time in between two consec- 
utive forward steps whereas r is the dwell time in be- 
tween two consecutive backward steps. Similarly, is 

the dwell time in between two consecutive steps of which 
the first is forward and the second is backward whereas 
the opposite if true in case of r |_. 

We denote the probability density functions for the 
conditional dwell times by the symbols 3±± (t). The inte- 
grated probabilities obtained from these probability den- 
sities are given by 



V±±(t)= f3 ±± {t')dt' 
Jo 



(31) 



Obviously, lim t ^ 00 V±±(t) = f*S±±(f)df = 1. More- 
over, we introduce the "pairwise splitting" probabilities 

Tl±± where n++ and 11^ represent the probability that 

a forward step is followed by a forward step or a backward 

step, respectively. Similarly, II_ + and II denote the 

probability that a backward step is followed by a forward 
step or a backward step, respectively. For the analysis 
of the experimental data and comparison with theoreti- 
cal predictions, it is sometimes more convenient to divide 
the dwell times into two groups depending on the direc- 
tion of the following step. We use the symbols S?j_ (t) and 
SI (t) to denote the probability density functions for the 
dwell times before a forward (+) and a backward (— ) 
steps, respectively. Obviously, S±(i) are given by @ 



_(t) = n ++ s ++ (t) + n_ + s + _(t) 
it) = n+_s_+(t) + n__s__(t) 



(32) 



In this section we derive analytical expressions for 
H±±, H±±(i), and hence, 3±(t) following the procedure 
adopted in ref . [|, [Toj] . 

In our model, immediately after a forward step the 
KIF1A motor can be found in either of the two states 
{Qn(t)jQi{t) = S(t) i.e. strongly bound state, and 
Q2{t) = W{t) i.e. weakly bound state) whereas it can 
exist only in state 2 immediately after a backward step. 
The first escape problem for our model is governed by a 
reduced Master equation (as described in refs. [|| and 
fl7j). with the transition matrix 



T = 



-UJ h U! s 

u h -{ojf + u s + 2u) b ) 



The eigen values of T are 

— Ul ± y/ (UJh - - Ul s - 2uj b ) 2 + 4bJ h UJ s 



Ai 



(33) 



(34) 



where, u\ = ui^ + ujf + u> s + 2uj b and define u 2 — AiA 2 = 
id f ujh + 2u)bU>h- The initial conditions, describing the dis- 
tribution of states just after a ± (forward and backward) 



step, are given by 



;(o) 



U) b + U)f 



andQ-(0)=^ 2 (35) 



where /i = 1, 2, Qi = S and Q 2 = W as described above. 

Using the expressions for the probability currents as- 
sociated with the allowed transitions in our model, we 
get 

n±+V±+(t) = f dt'{uj b Q 2 (t')+uj f Q 2 {t')} 
Jo 

n±-P±-(t) = / dt'{u b Q 2 (t')} (36) 
Jo 

Following Linden and Wallin [i| , by using the Ansatz 



V ++ (t)=V—(t) = l + J2^ 



,A„t 



(37) 



together with eqs. (f3"5" j) - (j3l)j) to compute d t H±+V±+(t) 
and dtH±-V±-(t) we obtain the following systems of lin- 
ear equations: 



(38) 





fet 


f +- 

6 1 


fc l 




Ai A 2 ) 1 




fc 2 


fc 2 




-1 


-1 




1 


-1 




2 
<*>b 




\-Ulf 




n++ n + _ 




Uf) n 


- + 





Solving this and using H±±(i) = dtP±±(t) we obtain 
the following distributions of the conditional dwell times: 



Hence 



= u 2 


( e A lt _ 


e A a tJ 


Ai - 




= u 2 


(e>it _ 


J + 


Ai- 




= u 2 


(gAit _ 


e A 2 t) 


Ai- 




= u 2 


( e Ait _ 




Ai- 




I n 2 


- n 1+ 

- n 2 + 


)-( 



=Ait _ \ c A 2 n 



n_+ (A1-A2) 



n 



(Ai - A 2 ) 





UJ b (ujf + UJ b ) 



(39) 



(40) 



To derive the pairwise splitting probabilities we first 
solve for the Hfj,± and then weight them according to the 
initial conditions as of eq. (j3"5")) : 

n ±+ = ^n, + Q( ± )(o),D ± . =^n^Qf)(o) (4i) 

H fj, 

We finally get the following splitting probabilities 

u) h (ujf +0J b ) 



n++ = n_+ = 
n = n+_ = 



u) h {u) f + 2ui b ) 
u 2 



(42) 
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FIG. 3: The distribution E5_(t) is plotted for few different 
values of 

Thus, probability of a forward step is [wh{ujf + 
W5)]/[w),(w/ + 2wh)] irrespective of the direction of the 
preceeding step. Similarly, the probability of a backward 
step is ((jJhk>b)/[wh{ U) f + 2u)} ) )] irrespective of the direction 
of the preceeding step. Both these are consistent with the 
kinetic pathways shown in fig. Q] as well as with the fact 

that LT ++ + LT + _ = 1 and II + = 1. 

Substituting the expressions PS|) and into (j3"2")l 
we get the analytical expressions for E±(t). This distri- 
bution is plotted in fig. [3] for a few different values of 
LOb- The most probable dwell time before a forward step 
decreases with increasing w&. 

V. RANDOMNESS PARAMETER 



FIG. 4: The randomness parameter r, defined by eq. (|46[) . is 
plotted against ATP concentration for a few different values 
of a;;,. The inset shows the same for few different values of 
w/. 
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Using eq. (|30|) we obtain 

(r 2 + ri)(r 2 - ri) 



< t >= 4a 



and 



< t z >= 16a 



(rir 2 ) 2 

{r 2 - ri)(r% + nr 2 + r|) 
(rir 2 ) 3 



(43) 



(44) 



where r 1; r 2 are given by eqs. (fM|) . (f2l)j) and a is given 
by 



2w 2 + u b ujf + u b uj s - uj b uj h - ujfUj h - uj b r 



ro 



(45) 



ro is given by eq. (j!?Tj) . Using equations (P), (|4^)) and flU]) 
we obtain randomness parameter as follows 

r _ r\r 2 (r\ + rir 2 + rf) - a(r2 - ri)(r 2 + ri) 2 
a(r- 2 - ri)(ri + r 2 ) 2 

Randomness parameter (see eq. (|4"?J]) ) is plotted in fig. |U 
against ATP concentration for a few different values of 
uib and same is plotted in the inset for a few different 



FIG. 5: The randomness parameter r+, defined by eq. (|47[) . 
is plotted against ATP concentration for a few different val- 
ues of LJb- The inset depicts the dependence of randomness 
parameter r*_, defined by eq. (|47p . against ATP concentration 
for the same set of tot values. 



values of ojf. With increasing concentration of ATP, r 
decreases and finally sturates near to unity. At low ATP 
concentrations r is greater than 1 which may be the effect 
of multi-exponcntiality in dwell time distribution. It also 
depicts that increasing concentration of ATP reduces the 
fluctuations in dwell time. 

For the conditional dwell times, the randomness pa- 
rameters are defined by 



< (4) 2 > - < t*± > 2 
<#. > 2 



(47) 



Using eq. (|32)l we obtain 

u 2 (Ai + A 2 ) 



< tl >- 



-2n, 



n+W 



(AiA 2 ) 2 II + _(w b + cj/)AiA 2 A X A 2 

(48) 
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FIG. 6: The randomness parameter r+, defined by eq. (|47|l . 
is plotted against ATP concentration for a few different val- 
ues of ujf. The inset depicts the dependence of randomness 
parameter r*_, defined by eq. (|47p . against ATP concentration 
for the same set of uif values. 



and 

W , M 2 S Arr u 2 (\l + X^ 2 + \ 2 2 ) 2 Wb (A!+A 2 ) 
< (t + ) > — 4H++- 



-2n 



++ 



(A1A2) 3 (AiA 2 )2 
w b 2 (Ai + A 2 ) 



n+_(w & +o; / )(AiA 2 ) 2 



(49) 



w 2 (Ai+A 2 ) n + _ (uib w 6 



where , II h ,II__, and 11^ are given by eq. (fl2"1) . 

Ai and A 2 are given by equation (1341) . Similarly, 



< t*_ >= -2n__ 

and 

1* \2 



n_+AiA 2 AiA 2 
(50) 



< (t* + y > = 4n__ 



-2n + _ 



(AiA 2 ) 2 

u 2 (Xl + X1X2 + X 2 ) , 2w b (Ai + A 2 ) 



(AiA 2 )3 
((j b +o;/)(Ai + A 2 ) 
n_+(AiA 2 ) 2 



(AiA 2 ) 



(51) 



The conditional randomness paremeters (eq. (|47p ) are 
plotted in the figs. [5] and [5] against ATP concentration 



for different values of u>b and w/ respectively. The non 
monotonic variation of conditional randomness parame- 
ters with ATP concentration changes to a monotonic de- 
crease (see fig. E]) when the magnitude of LOb is sufficiently 
high. Variation of the randomness parameter indicates 
the changes in the number of rate-limiting steps. Any 
value of the randomness parameter higher than unity 
may, at first sight, appear counter-intuitive. But, this 
is quite common for systems with branched mechano- 
chcmical kinetics. 



VI. SUMMARY AND CONCLUSION 

Theoretical calculation of the dwell time distribution 
of two-headed conventional kinesin motors has been re- 
ported earlier [12[ . In this paper we have derived an exact 
analytical expression for the distribution of dwell times 
of single-headed kinesin motors KIF1A at each binding 
site during a processive run on a microtubule. We have 
used the NOSC model [| for single-headed KIF1A mo- 
tors to derive our results. Since both forward and back- 
ward steps of this motor are possible, we have also de- 
fined conditional dwell times and calculated their distri- 
butions analytically. The experimentally measured dwell 
time distributions for some of the other processive mo- 
tors, like conventional kinesin (l3| . myosin- V [l4j . dynein 
[l5l |. ribosome !16], etc., have been reported in the lit- 
erature. To the best of our knowledge, the dwell time 
distribution of KIF1A has not been reported so far; we 
hope our theoretical prediction will stimulate experimen- 
tal investigations. 
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